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ABSTRACT 

We study internal extinction of late-type galaxies in the Sloan Digital Sky Survey. We find that 
the degree of internal extinction depends on both the concentration index c and Xg-band absolute 
magnitude Mk- We give simple fitting functions for internal extinction. In particular, we present 
analytic formulae giving the extinction-corrected magnitudes from the observed optical parameters. 
For example, the extinction-corrected r-band absolute magnitude can be obtained by Mr^ = —20.77-1- 
(-1 + ^\ + M\{M^^ohs + 20.77 + 4.93A)) /2A, where A = 0.236{1.35(c - 2.48)^ - 1.14} log(a/5), c = 
-Rgo/^so is the the concentration index, and ajh is the isophotal axis ratio of the 25 mag/arcsec^ 
isophote in the z-band. The Icr error in M^fi is 0.211og(a/6). The late-type galaxies with very 
different inclinations are found to trace almost the same sequence in the (u — r)-Mr diagram when 
our prescriptions for extinction correction are applied. We also find that (u — r) color can be a third 
independent parameter that determines the degree of internal extinction. 

Subject headings: galaxies: general-galaxies:fundamental parameters-galaxies: ISM -galaxies: spiral- 
(ISM:) dust, extinction 



1. INTRODUCTION 

Dust in spiral galaxies causes internal extinction. For 
example, due to dust, edge-on spiral galaxies in general 
look fainter than face-on spirals with the same intrin- 
sic luminosity in short wavelength optical bands, in par- 
ticular (see Fig. 12 of Choi, Park, & Vogeley 2007). 
The study of internal extinction is important for deter- 
mination of distances and absolute magnitude of galax- 
ies. Giovanelli et al. (1995) demonstrated that inade- 
quate treatment of internal extinction has strong effects 
on the distances estimated using the TuUy-Fisher relation 
(TuUy & Fisher 1977). Internal extinction also causes bi- 
ased measurements of galaxy star formation rates (see, 
for example. Bell & Kennicutt 2001; Sullivan et al. 2000) 
and affects determination of galaxy luminosity function 
(Shao et al. 2007). Earlier studies of internal extinction 
include Giovanelh et al. (1994, 1995), TuUy et al. (1998), 
and Masters, Giovanelli, & Haynes (2003). Recent works 
that contain relevant discussions on the topic are Rocha 
et al. (2008), Unterborn & Ryden (2008), and Mailer et 
al. (2008). 

The amount of obscuration by dust is larger when the 
observing wavelength is shorter. Therefore, the effect 
of inclination is more pronounced in short-wavelength 
bands, such as u or r band, and it may be much smaller 
in Ks (2.17/im) band. As a consequence, when the view- 
ing angle changes for a given galaxy, u- or r-band magni- 
tude changes while .ftTs-band magnitude does not change 
much. Therefore, u — Kg or r — Kg color tends to be 
larger when the galaxy is viewed more edge-on. 

Internal extinction depends on many factors. Perhaps, 
the most important factor is the amount of dust. How 
dust is distributed can also be an important factor. The 
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study of extinction versus inclination will ultimately re- 
veal how dust is distributed and how much dust is con- 
tained in spiral galaxies. Then, what determines the 
amount and distribution of dust in spiral galaxies? 

Earlier works have discussed dependence of internal 
extinction on luminosity. Giovanelli et al. (1995) found 
that the amount of internal extinction in / band depends 
on the galaxy luminosity. TuUy et al. (1998) also found a 
strong luminosity dependence using a magnitude-limited 
samples drawn from the Ursa Major and Pisces Clus- 
ters. Masters, Giovanelli, & Haynes (2003) studied in- 
ternal extinction in spiral galaxies in the near infrared 
and also found a luminosity dependence. On the other 
hand, there are suggestions that internal extinction de- 
pends on galaxy type (de Vaucouleurs et al. 1991; Han 
1992). 

In this paper, we study the internal extinction in SDSS 
late-type galaxies. We investigate the dependence of the 
inclination effects on luminosity, concentration index c, 
and u — r color. In §2, we describe the data set used in 
this paper. In §3, we study dependence of internal extinc- 
tion on the concentration index c and i^s-band luminos- 
ity separately. In §4, we measure dependence of internal 
extinction on the concentration index c and Kg-hand lu- 
minosity simultaneously. In §5, we derive dependence of 
inclination effects on r-band luminosity. The result in §5 
is useful when X-band magnitude is not available. In §6, 
we present dependence on u — r color. We give discussion 
in §7 and conclusion in §8. 

2. DATA AND METHOD 

In this study, we investigate how u- and r-band magni- 
tudes behave as the inclination angle changes. The phys- 
ical parameters we consider are r-band absolute magni- 
tude Mr, u — r color, the axis ratio a/5, and the concen- 
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tration index c. 

The primary data set we use is a subset of volume- 
limited SDSS (DR5) galaxy sample (the data set Dl in 
Choi et al. 2007). The redshift of the galaxies is be- 
tween 0.0250 and 0.04374 and the minimum r-band ab- 
solute magnitude, Mr, is —18.0 + 51og/i, where h is the 
Hubble constant divided by 100 km s^^ Mpc^^. In this 
paper, we assume the Hubble constant is 75km/sec/Mpc. 
Therefore, the absolute magnitude in this paper can be 
transformed to the /i-dependent form by 

Mi^ = M^=°-^5 + 51og(V0.75), (1) 

where M^^^-''^ is the absolute magnitude we use in this 
paper. The rest-frame absolute magnitudes of galaxies 
are computed in fixed bandpasses, shifted to z = 0.1, us- 
ing Galactic reddening corrections (Schlegel, Finkbeiner, 
& Davis 1998) and iif-corrections as described by Blan- 
ton et al. (2003). Therefore, all galaxies at 2; = 0.1 
have a _ft'-correction of — 2.51og(l -I- 0.1), independent of 
their spectral energy distribution. We then apply the 
mean luminosity evolution correction given by Tegmark 
etal. {2QQA),E{z) = 1.6(^-0.1). The comoving distance 
limits of our volume-limited sample are 74.6 and 129.8 
ft,~'^Mpc if we adopt a flat ACDM cosmology with den- 
sity parameters 0^ = 0.27 and = 0.73. The galaxies 
in the volume-limited sample are divided into early (E 
and SO) and late (S and Irr) morphological types based 
on the location of galaxies in the u — r color, g — i color 
gradient, and concentration index space (Park & Choi 
2005). 

There are 14,032 late- type galaxies in the data set and, 
among them, ~8,700 galaxies have matching data in the 
2 Micron All-Sky Survey (2MASS) Extended Source Cat- 
alog. If the distance between the center of a SDSS galaxy 
and that of a 2MASS galaxy is less than the semimajor 
axis in the SDSS i-band, we consider they are identical. 
When there are multiple matches, we simply discard the 
data. We also remove data when the absolute value of 
the color gradient parameter (A(g — i); see Park & Choi 
2005) is larger than 0.7, the u — r color is larger that 4 or 
less than 0, or the Hq line width is less than or larger 
than 200A. About 700 galaxies have been discarded from 
this procedure. 

2.1. Parameters 

1) Mr, Mu, and u — r: We use r-band absolute Pet- 
rosian magnitudes and u — r model color from the SDSS 
data. We obtain u-band absolute magnitude M„ from 
the relation 

Mu = Mr + {u-r). (2) 

Since the uncertainties in the extinction-corrections are 

large, we ignore the differences between the Petrosian 
u-band absolute magnitude and that derived from Equa- 
tion (2). 

2) c: The concentration index is defined by Rqq/R^q 
where J?5o and Rqq are the semimajor axis lengths of 

ellipses containing 50% and 90% of the Petrosian flux 
in the SDSS i-band image, respectively. It is corrected 
for the seeing effects by using the method described in 
Park & Choi (2005). It is basically a numerical inverse 
mapping from the image convolved with the PSF to the 
intrinsic image having a Sersic profile and a fixed incli- 
nation. The concentration index used here is the inverse 
of the one used by Park & Choi. 
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Fig. 1. — Distribution of the galaxies in our volume-limited sam- 
ple in the ii's-band absolute magnitude Mk and the concentration 
index c = Rqo/Rtm space divided into 25 X 25 cells. 

3) a/b: The isophotal axis ratio a/& is from the SDSS 
i-band image. Here, a is the major axis length and b 
the minor axis length, corrected for the seeing effects. 
We choose the isophotal axis ratio because the isophotal 
position angles correspond most accurately to the true 
orientation of the major axis, which is true for barred 
galaxies in particular. Here we assume that outside the 
central region, the disk of a late-type galaxy can be ap- 
proximated as a circular disk, which thus appears as an 
ellipse in projection. If the disk of late-type galaxies has 
an intrinsic non-circularity, our estimation of the a/b ra- 
tio will have some error due to our assumption. 

4) Mk'- We use i^^-band magnitude (more pre- 
cisely, magnitude within the 20th mag arcsec^^ elliptical 
isophote set in Kg band) from the 2 Micron All-Sky Sur- 
vey (2MASS) data. We also apply i^sT-corrections and 
Galactic extinction corrections for the 2MASS data as 
described by Masters, Giovanelli, & Haynes (2003). As 
in the SDSS case, the if-corrections are made for a fixed 
redshift of 2 = 0.1. Most matched galaxies in our sample 
have Mk brighter than —20. 

2.2. Method 

We adopt the popular parameterization of the effects 
of internal extinction: 

= 7Alogio(a/&), (3) 

where A\ is the amount of extinction. In this parame- 
terization. 7a is the slope of a scatter plot drawn on the 
logio(a/fe) - A\ plane. We will see in section 4 that this is 
actually a very reasonable model in the case of the SDSS 
galaxies even though some recent studies reported that 
better models can be found (Masters et al. 2003; Rocha 
et al. 2008; Unterborn & Ryden 2008). 

Our goal is to find the values of 7^ and 7„. To find ^r 
(or 7^), we first plot r-Kg (or u-Kg) against \ogiQ{a/b). 
Then, the slopes of the scatter plot are 

jr-K^Jr -JK, and 

lu-K^-yu-jK- (4) 

When no strong extinction is present in Kg band, we can 
write 

7r-K~lr, and 

7«-if~7«- (5) 

In this paper, we assume that Xs-band magnitude is al- 
most free of inclination effects. Masters et al. (2003) 
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analyzed galaxies in the 2MASS Extended S ource Cata- 
log and concluded that the internal extinction is indeed 
small in Kg band. Their results are consistent with the 
earlier result that 7^ ~ 0.22 (TuUy et al. 1998). 

As shown by earlier studies, 7a depends on luminosity 
(Giovanelh et al. 1995; Tully et al. 1998) and/or galaxy 
type (de Vaucouleurs et al. 1991; Han 1992). In this 
study, our primary concern is the dependence of 7^ on 
i^Tg-band luminosity and the concentration index c. In 
Figure 1, we plot the galaxies on the c — Mk plane. On 
the plot, we show 25 x 25 cells. The coordinates of the 
cell centers, {cm,Kn), are given by 



c,n = 1.68 + 0.06m, m= 0, 1,..., 24, and 



-23.4 + 0.15n, 



n= 0, 1,..., 24. 



(6) 
(7) 



Note that C2 = 1.8, C22 = 3.0, K2 = -23.1, and ^^22 = 
-20.1. 

When we investigate the dependence of 7^ on Mk and 
c, we may simply use galaxies in each cell in Figure 1. 
However, some cells in Figure 1 are not sufficiently pop- 
ulated. Therefore, for smooth results, we use galaxies 
in 5 X 5 cells for scatter plots. More precisely, we use 
galaxies with Cm — 2. 5 (Ac) < c < Cm + 2. 5 (Ac) and 
Kn-2.5{AK) <Mk <K„ + 2.5iAK),wheie Ac = 0.06 
and AK = 0.15. We allow for a similar overlap when we 
study Mk- or c-dcpendence separately. 

To find the slopes, we divide log(a/6) axis into 20 bins 
between and 1. Then we find a representative value for 
each bin. We try two methods to obtain the representa- 
tive value: 

1) the median 

2) the most probable value. 

To fine the most probable value in each log(a/6) bin, we 
use a smoothing function 

,A_/exp(-(2/-j/i)V2cT^) if\x-Xi\<Om, . 
nx,y)-<^Q otherwise, 

where x = logio(a/6), Xi {yi) is the x (y) value of the i*^ 
galaxy, and ay is the standard deviation. After applying 
the smoothing function, we find the maximum value of 
the smoothed distribution in each bin. After finding rep- 
resentative values in bins of log(a/6), we perform the lin- 
ear fit to the representative values between log(a/6) = 
and 0.5. 

3. DEPENDENCE ON C OR Mk 
3.1. Dependence on c 
In this subsection, we consider only c-dependence. 



We plot r - K, 



Kg, and u — r colors against 



logiQ{a/b) in Figure 2. As we mentioned earlier, when 
we draw the scatter plot for c = Cm , we use galaxies with 
Cm - 2.5(Ac) < c < Cm + 2.5(Ac), where Ac = 0.06. 
Therefore, the plot for c = 1.80 contains galaxies with 
1.65 < c < 1.95, for example. 

We can see that the slope, hence 7^, for c = 1.80 is 
smaller than those for other c values. Note that the slope 
for c = 1.80 is very close to zero for u — r color. 

Figure 3 shows dependence of the slope on c. All 3 col- 
ors show a common feature: the slope peaks near c ~ 2.5. 
This means that internal extinction is maximum for in- 
termediate late-type galaxies, and is smaller for early 
and late late-type galaxies. On the other hand, u — r 



and u — Ks show stronger dependence on c than r — Ks, 
telling that internal extinction is higher in shorter wave- 
length bands. The quadratic equations on the plots are 
the fitting functions. The solid curves are for the most 
probable values and the dotted curves for the median. 
The RMS scatters of the measured slope from the fitting 
functions (for the most probable values) are 0.071, 0.170, 
and 0.154 for r — Ks,u — Kg, and u — r colors, respec- 
tively. We list the quadratic fits, 7(c) 's, in Table 1 and 
2. 

3.2. Dependence on Mk 

In this subsection, we consider M^-dependence of in- 
ternal extinction. We plot r — Kg, u — Kg, and u — r 
colors against \ogiQ{a/b) in Figure 4. When we draw 
the scatter plot for Mk = Kn, we use galaxies with 
K,, - 2.5{AK) < Mk < K„ + 2.5(Air), where AK = 
0.15. For example, the plot for Mk = —23.1 contains 
galaxies with -23.475 < Mk < -22.725. 

We can clearly sec that the slope, hence 7^, for Mk = 
—20.1 is smaller than those for other Mk values. Note 
that the slope for Mk = —20.1 is very close to zero for 
all 3 colors. 

Figure 5 shows dependence of the slope on Mk- All 3 

colors have a common feature: the slope is higher when 
galaxies are brighter in Kg band. Due to lack of data 
points, the slope is not clear for Mk < —23.2. There- 
fore, one should be careful when using the fitting func- 
tions shown on the plots for galaxies with Mk < —23.2. 
The RMS scatters of the measured slope from the fitting 
functions (for the most probable values) are 0.071, 0.110, 
and 0.075 for r — Kg, u — Kg, and u — r colors, respectively. 
We list the quadratic fits, 7(Mx)'s, in Table 1 and 2. 

4. DEPENDENCE ON C AND Mk 

We now study dependence of r — Kg, u — Kg, and u — r 
colors on the concentration index and Kg-h&nA absolute 
magnitude. In Figure 6 we plot u — r against \ogiQ{a/h) 
as a function of c and Mk- The diamonds are the most 
probable values in bins of log(a / h) . The lines are the least 
square fits to the most probable values. It can be seen 
that the reddening is fit well by our extinction model 
linear in logiQ(a/6) (i.e. Equ. 3), and that the slope 
depends on both c and Mk- Plots for r — Kg and u — Kg 
show behaviors similar to the u—r case. We also obtained 
the fits to the median, which are qualitatively similar to 
the most probable case. 

In the top panels of Figure 7 we present contour plots 
of the measured slopes oir — Kg color, ^r-K- The upper- 
left panel is for 7r-x based on the most probable values. 
The function 



7^_k(c,Mk) = l-^2-ir-K{c)lr-K{MK), 



(9) 



shown in the upper-middle panel, fits well the slopes 
based on the most probable values (see Table 1). Since 
we do not have enough bright galaxies in Kg band, we 
do not have a reliable fitting formula for Mk < —23.25. 
Although our fitting formula suggests that the slope de- 
clines when Mk becomes less than ~ —22.9, the true 
behavior of the slope may be different. Therefore, in- 
stead of using the fitting formula in Equation (9), one 
may use 



7^_k(c,Mk) = 1.02 X 1.147^_k(c) 



(10) 
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c=1.80 c = 2.10 c = 2.40 c = 2.70 c = 3.00 
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Fig. 2. — Dependence of colors on the concentration index c (= Rgo/Rso)- When c is small, the colors only weaJjly depend on the axial 
ratio a/6. The axial ratio is derived from the SDSS i-band image. 




Fig. 3. — Dependence of the slope of the extinction law on c {= Rgo /R50). Solid lines are the slopes of the linear fit to the most probable 
values in log(a/6) bins between log(a/b) = and 0.5. Dotted lines are the slopes of the linear fit to the median. {Left panel) 7, — x- {Middle 
panel) 'Ju-K- {Right panel) 'yu-r- 



for Mk < —22.9. The upper-right panel of Figure 7 
shows the difference between the measured slope and the 
fitting function. The RMS value of the difference is 0.174 
while the peak slope is about 1.4. See Table 2 for a fitting 

function based on the median. 

Similarly, wo present contours of 7u_/i: and ju-r in the 
middle and bottom panels of Figure 7, respectively. The 
fitting functions for the most probable values are 

^u-k{c,Mk) = 0AS-1u-k{c)iu~k{Mk), (H) 
7„_,(c, Mk) = 0.947„_,(c)7„_,(Mk) (12) 
(see Table 1). Again the fitting functions are uncer- 
tain for Mk < -23.2. The RMS differences, shown in 
the middle-right and lower-right panels of Figures 7, are 
0.354 and 0.234, respectively. See Table 2 for fitting func- 
tions based on the median. 

5. DEPENDENCE ON r MAGNITUDE 



When i^s-magnitude is not available, we cannot use 
the fitting functions derived in the previous section. 
In this section we present a method that utilizes r- 

magnitude, instead of ii's-magnitude. We consider 7's 
derived from the most probable values. 

5.1. Dependence on M^.o 

We obtain intrinsic (or face-on) r-band absolute mag- 
nitude, Mj-fl, from 

Mr,o = Mr,obs - 7r{c, Mk) logio(a/6), (13) 

where M,..ohs is the r-band absolute magnitude before 
inclination corrections and 7r-(c, M^) ~ 7r-K(c, M/f) is 
given in Table 1. Note that 7,._if(c, M^) > 0. The 
RMS differences between the measured slopes and the 
fitting functions given in the last column of Table 1, are 
0.208, 0.390, and 0.233 for r - Ks,u - K^, and u - r 
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Fig. 4. — Dependence of colors on Mk- Faint galaxies in Ks band exhibit shallower slopes in all 3 colors, while bright galaxies show 
steep slopes. 




Fig. 5. — Dependence of the slope of the extinction law on the -ftTs-band absolute magnitude Mk- Solid lines are the slopes of the linear 
fit to the most probable values in log(a/6) bins between log(a/6) = and 0.5. Dotted lines are the slopes of the linear fit to the median. 
{Left panel) 7, — k- {Middle panel) 7u_k- {Right panel) 'yu-r- 



colors, respectively. After obtaining Af,.,o we investigate 
how internal extinction depends on the intrinsic r-band 
luminosity. 

Scatter plots in Figure 8 clearly shows that the slopes 
of the extinction in r — Kg, u — Kg, and u — r colors de- 
pend on the r-band absolute magnitude M^fi. The slopes 
of the scatter plots are very small when M^fi = —18.80, 
which means there is very little extinction in galaxies 
much fainter than the M* galaxies. (But note that 
late-type galaxies tend to be irregulars at fainter mag- 
nitudes and that irregulars have patchy distribution of 
dust which can cause a weaker dependence of extinction 
on inclination.) 

Figure 10 show the behavior of the slopes on c — Mr,o 
plane. It is interesting that there is a well-defined peak 
near (c, Mr ft) ^ (2.5, —20.8). In the r-band, the internal 
extinction is maximum at the absolute magnitude very 



close to the characteristic magnitude (see Table 2 of 
Choi et al. (2007) for measurements for the late- 
type SDSS galaxies. Note that h = 0.75 is used in the 
present work.) The contour spacing for -yr-K (left panel) 
is 0.2 and the value of 7^-^- at the center of the peak is 
greater than ^ 1.6. We can also see a similar peak for 
^u-K (middle panel). However, the location of the peak 
for 7„_r (right panel) is somewhat different. 

Figure 9 shows the dependence of the slope on Mr,o- 
The fitting functions are also shown on the plots. 

5.2. Obtaining Mr,o from Mr,obs and c 

We cannot observe the intrinsic r-band absolute mag- 
nitude Mj.fi directly. But we can obtain M^fi by solving 
the following equation: 

Mrfi + 7r-K{c, Mrfi) logio(a/6) = Mr^obs, (14) 
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Fig. 6. — Scatter plots for u — r. Diamonds denote the most probable values in log(a/6) bins. Lines are the results of the least square 
fit to the diamonds. Only points with logi()(a/6) < 0.5 are used for fitting. 
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where Mr.obs, c, and a/b are observed quantities. When 
we use a quadratic approximation as shown in Figure 9, 
then Eq. 14 becomes 

Mr,o + A[{Mr,o + 20.77)2 - 1.10/0.223] = Mr,obs, (15) 
where 

A = 1.06 X 0.223 X [1.35(c- 2.48)2 -1.14] logio(a/6) < 0. 

(16) 

Equivalently. wc have 

(Mr o+20.77)+A[(M^ 0+20.77)2-4.93] = (M^ o&s+20.77). 

(17) 

The solution is 

, -1 + V^l + 4A{Mr,obs + 20.77 + 4.93A) 

J^^r,0 = —-i^-l l-\ ^ , 

which gives the correct answer when A goes to zero. This 
equation is valid for [{Mrfl + 20. 77)^ - 4.93] < and 
A < 0. When [(M^^o + 20.77)^ - 4.93] > or A > 0, we 
simply have M^^o = Mr^obs- 

Figure 11 shows that extinction corrections using Mr 
and c give results quite close to those using Mk and c. 
In the left panel, we obtain the intrinsic r-band absolute 
magnitude using two methods. For the x-axis, we make 
inclination corrections using 

7r ~ 7r-K = 7r-if (C, Mk), (19) 

where jr-KicjMx) is given in Table 1. For y-axis, we 
make inclination corrections using Eq. (18). The left 
panel shows that the two methods give a good agreement. 
The middle and right panels arc obtained similarly. The 
la error in Mrfl derived from Mr^obs and c is given by 
the error in 7r-A:(c,A'/r g) times log(a/6) (see Eq. 14) 
or 0.2081og(a/6) over the range from M^o = —18.5 to 
-22.5. 

5.3. Obtaining M^^q and {u — r)o 

We can estimate M„^o from Mj-^obs and c. First, we 
need to obtain M„^o as described in the previous subsec- 
tion. Mtj.o is 

M„,o = M^.ob,, - 1u-k{c, Mr,o) logio(a/6), (20) 

where Ju-k{c, Afr,o) is given in Table 1. 

Wc can also estimate {u — r)o from a/b and c. First, 
we need to obtain Mrfl. Then, (u — r)o is 

{u - r)o = {u- r)a/b - 7«-r(c, Mr^) logio(a/6), (21) 

where ju-r{c, Mr^) is given in Table 1. 

6. DEPENDENCE ON t/ - i? COLOR 

In this paper, we have mainly considered dependence 
of 7a on K and c. There may be more parameters that 
determine 7^. In this section, we show that u — r color 
can be a third independent parameter. 

6.1. Slope vs. {u — r)o 
We obtain the intrinsic color, {u — r)o, as follows: 
{u - r)o = {u- r)a/b - 7u-r(c, Mrfi) logio(o/&), (22) 
or, 

{u - r)o = {u- r)a/b - 7«-r(c, Mk) logio(a/&), (23) 



where "fu-r{c, Mr_o) and "/u-r{c, Mk) are given in Table 
1. In this section, we use Eq. (23). 

Figure 12 shows that -fr-K strongly depends on (u — 
r)o. The Figure shows that the slope for r — Kg is vir- 
tually zero when (u — r)o is less than ~ 1.0. The slope 
increases as the color increases. It seems that the slope 
reaches a maximum at {u — r)o ^ 2.5. The behavior of 
the slope is uncertain for {u — r)o > 2.5. 

6.2. Dependence of {u — r)o on c or Mk 

In this section, we investigate the possibility of (m — r)o 
as a third parameter. The third parameter should be in- 
dependent of the first 2 parameters (i.e. c and Mk in 
our case). However, Figure 6 shows that {u — r)o, which 
corresponds to the y intercepts of the linear fits, depends 
on the both c and Mk- Therefore, {u — r)o is not a com- 
pletely independent parameter, which means we cannot 
tell whether the change in slope observed in Figure 12 
is entirely due to {u — r)o or just a different realization 
of c or Mk dependency through the correlation between 
them and {u — r)o. 

However, a careful look at Figure 6 reveals that {u — r)o 
varies between ~ 1.5 and ~ 2.0. One exception is the 
{u — r)o (i.e. the y intercept) in the upper-right panel, 
which is as large as ^ 2.4. From this observation, we 
can conclude that the intrinsic scatter in (u — r)o itself is 
larger than the systematic change in (u — r)() due to c or 
Mk change. Therefore, if we properly limit the ranges of 
c and Mk and draw a plot similar to Figure 12, then we 
can tell whether or not ^/r-K really depends on {u — r)o. 

For this purpose, we only consider galaxies in the fol- 
lowing ranges: 

group 1: -22.0 < Mk < -21.4 and 1.8 < c < 2.1 
group 2: -22.0 < Mk < -21.4 and 2.1 < c < 2.4 
group 3: -21.4 < Mk < -20.8 and 1.8 < c < 2.1 
group 4: -21.4 < Mk < -20.8 and 2.1 < c < 2.4 

Then, we study the dependence of jr-K on (u — r)o for 
each group. 

6.3. {u — r)o color as a third parameter 

The discussion above indicates that the systematic 

change in {u — r)o in each group due to changes in c 
and Mk should be marginal. Figure 13 shows the possi- 
bility that {u — r)o can be a 3rd parameter: For given c 
and Mk, the slope 7^ shows dependence on (m — r)o. The 
general trend is that the slope is steeper, when (u — r)o 
is larger. Note, however, that galaxies in group 4 do not 
show strong dependence on {u — r)o, while those in group 
3 show a very strong dependence. 

7. DISCUSSION 

Qualitatively speaking, our results are consistent with 
earlier claims that the extinction in spirals depends on 
luminosity (Giovanelli et al. 1995; Tully et al. 1998). 

Tully et al. (1998) derived an extinction law in R band. 
They found that galaxies with Mr reaching —21 mag 
have "fn ^ 1.16 and the amplitude of extinction drops 
rapidly as luminosity decreases. Similarly, our Figure 9 
shows that "fr-K ~ 1-1 for Mr ~ —21 and it declines 
rapidly with decreasing luminosity. However, there are 
also differences. For example, "fn in Tully et al. (1998) 
reaches as large as ~ 1.4 for galaxies brighter than -21 
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y Fitting Ftn Difference 




Fig. 7. — The measured slopes of the extinction law (left panels), fitting functions (middle panels), and the difference between the 
measured slopes and the fitting functions (right panels). {Upper panels) r-K^ color. {Middle panels) u-K^ color. {Lower panels) u — r 
color. The horizontal axis of each contour plot is the concentration index c (= R90/R50) and the vertical axis is the Kg-hand absolute 
magnitude Ma- The contours are drawn at intervals of 0.2 (upper panels), 0.4 (middle panels), and 0.2 (lower panels), respectively. 



Internal Extinction 



9 



mag in the R band. We do not observe such a rise of 7^ in 
our results. However, since the definition of r-magnitudc 
is different and 7^ depends on other parameters, such as 
c and u — r, it is very difficult to understand the origin 
of this discrepancy. 

Shao et al. (2007) found 7^ ~ 1.37 when they used the 
r-band axis ratio, which is substantially larger than our 
values. The difference may stem from the difference in 
fitting methods. We studied how r — Kg color changes 
as a/b changes. Therefore the value of 7 in our study is 
actually 7r-if, not 7^. The value of 71. will be given by 
7r ^ Ir-K + IK ^ 1-15 + 0.26 ^ 1.41, whcrc 1.15 is the 
average value of jr-K near the maximum (see Figure 5) 
and 0.26 is the estimated (Masters et al. 2003). The 
definition of a/b is also different. Our a/b is the i-band 
isophotal axis ratio, while that of Shao et al. (2007) is 
the r-band axis ratio obtained from the best fit of the 
images of galaxies with an exponential profile convolved 
with the point spread function. 

On the other hand, there are several observations in / 
band. Giovanelli et al. (1994) found the extinction in I 
band of 7/ - 1.05 ± 0.08. Masters et al. (2003) foimd 
that 77 ~ 0.94. Therefore, our result of 7^. ~ 1.1 agrees 
with the common wisdom: there is more extinction in 
shorter passband. 

The concentration index, c, is related to morphology of 
late-type galaxies. Earlier studies show that, on average, 
early late-type galaxies have higher concentration (see, 
for example, Shimasaku et al. 2001; Strateva et al. 2001; 
Goto et al. 2003; Yamauchi et al. 2005). Therefore de- 
pendence of inclination effects on the concentration index 
implies that galaxy morphology is an important factor, 
which is consistent with earlier findings (see, for example, 
de Vaucouleurs et al. 1991; Han 1992). It is interesting 
that internal extinction is maximum when c ~ 2.5 which 
corresponds to Sb galaxies (see Shimasaku et al. 2001). 

Choi et al. (2007) have shown that the sequence of 
late-type galaxies in the color- magnitude diagram can be 
significantly affected by the internal extinction. When 
the sequences are compared between late types with 
b/a > 0.8 (nearly face-on) and b/a < 0.4 (nearly edge- 
on), a very large departure is observed for the sequence 
of the inclined late-type galaxies relative to that of the 
face-on late types as the inclined galaxies appear fainter 
and redder. The degree of departure was found to be 
maximum for intermediate luminosity late types with 
Mr ~ —20.5 -I- 51og/i = —21.1. No such extinction effect 
was found for early- type galaxies. It will be very inter- 
esting to see if our extinction-correction prescription re- 
stores the color-magnitude diagram of inclined late- type 
galaxies. 

In Figure 14, we plot our late-type galaxies in the color 
magnitude diagram. The left panel of Figure 14 shows 
the galaxies with a/b < 1.25 (b/a > 0.8) and the middle 
shows those with a/b > 2.5 {b/a < 0.4). We find that the 
group of highly inclined galaxies is significantly shifted 
with respect to the almost face-on ones toward fainter 
magnitude and redder colors. The result in Figure 12 in- 
dicates that the color of very blue galaxies is less affected 
by the internal extinction. However, Figure 12 does not 
tell much about color change of very red galaxies (i.e. 
redder than w — r ~ 2.5) due to lack of data. Note that 
Choi et al. (2007) observed that both very blue and very 
red galaxies are less effected by the internal extinction. 



The right panel of Figure 14 is the inclination-corrected 
color magnitude diagram. The inclination correction is 
done based on Mk- and c-dependence of Mr and u — r. 
To be more specific, we use Eqs. (13) and (21) to obtain 
Mrfi and {u — r)o. We do not draw the reddening vectors 
(caused by internal extinction) on the color-magnitude 
diagram bccasc, due to c-dcpendency, infinite number of 
reddening vectors arc possible for a given value of (Mr, 
u ~ r). If one wish to draw 'average' reddening vectors 
for a fixed value of log(a/6), one can do that by using 
7u_,.(Mr,o) and Jr{Mrfi) in Table 1 (or 2). Of course, 
to get the 7 values, we first need to obtain M^o from 
Eq. (18). 

8. CONCLUSION 
Our major conclusions arc as follows. 

1. We have shown that the slope 7^ depends on both 
Ks-hand absolute magnitude Mk and the concen- 
tration index c. The fitting functions for the rela- 
tions are given in Table 1 and 2. 

2. We have also shown that the slope 7a depends on 
both the concentration index c and r-band absolute 
magnitude Mrfl, where the subscript '0' denotes 
the value after inclination correction. The relations 
are also given in Table 1 and 2. 

3. We have derived analytic formulae giving the 
extinction-corrected magnitudes from the observed 
optical parameters (see, for example, Equation 
(18)). 

4. We have shown that (m — r)o can be a third param- 
eter that determines the slope 7a . 
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M^o=-21.80 M^o=-21.05 M^o=-20.30 M^o=-19.55 M,o=-18.80 




Fig. 8. — Dependence of colors on the intrinsic (i.e. face-on) r-band absolute magnitude M^-fi. We use ■yn ~ 7, — k ~ 1r-K{c, M^S) (see 
Table 1) to derive M^_o- When M^^o ^ — 19> internal extinction is very small. 




Fig. 9. — Dependence of 7 on the intrinsic (i.e. face-on) r-band absolute magnitude Mr.a. Solid lines arc the slopes of the linear fit to 
the most probable values in log(o/fe) bins between log(a/b) = and 0.5. Dotted lines arc the slopes of the linear fit to the median. The 
RMS scatters of the measured slope from the fitting functions (for the most probable values) are 0.064, 0.095, and 0.070 for r — Ks,u — Ks, 
and u — r colors, respectively. {Left panel) ^r-K- {Middle panel) 7u_if • {Right panel) -fu-r- 
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Ifr-K Ju-K ^u-r 




2.0 2.5 2.0 2.5 2.0 2.5 

c c c 



Fig. 10. — Dependence of the slope of the extinction law for the r — K, u — K, and u — r colors on the r-band absolute magnitude Mj.fi 
and the concentration c. The slope 7 is based on the most probable values. X-axis is the concentration index (= Rgo/Rso)- The contour 
intervals for r-Kg, u-Kg, and u — r are 0.2, 0.4, and 0.2, respectively. 




n,,„ from M,_„j. M.„ from M,_^. (u-r)„ from M,^, 



Fig. 11. — Comparisons between the inclination correction method using 7(0, Mk) and 7's in the shorter bands. The inclination correction 
using ■y{c,Mfc) is straight forward. However, the use of -f{c,Mrfi) requires one more step (see Eq. 18) to get Mrfi first. 



1.5 




(u-r)o 



Fig. 12. — Dependence of -yr-K on (m — r)o, where the sub- 
script '0' denotes the inclination corrected value. We use -fu—r = 
7u_r(c, M/f) (see Table 1) to derive the intrinsic color. 



Internal Extinction 



13 




Fig. 13. — {u — r)o as a third parameter. This figure shows 
that the slope depends not only on c and luminosity but also the 
intrinsic {u — r) color. Internal extinction is more significant when 
galaxies have larger intrinsic (u — r) color. 



a/b ^ 1.25 a/b > 2.5 all data 
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Fig. 14. — Distributions of latc-typc galaxies in the color magnitude diagram. Galaxies with a/6 < 1.25 {b/a > 0.8) arc shown in the 
left panel, and those with a/6 > 2.5 (6/a < 0.4) arc shown in the middle panel. Distribution of all galaxies after inclination corrections 
are drawn on the right panel. The median color-magnitude relations are also drawn. We find that the group of highly inclined galaxies is 
significantly shifted with respect to the almost face-on ones toward fainter magnitudes and redder colors. 
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TABLE 1 

Fitting Functions (for most probable values) 





Dependence on c 
7{c)^ 


Dependence on Mk 


Dependence on o 


& c 
7(c,Mk) 


Mrfi & c 
7(c,Mr,o) 


r-Ks 
u-Ks 
u — r 


-1.35(c-2.48)2+1.14 
-3.37(c-2.57)2+2.61 
-2.02(c-2.63)2+1.49 


-0.089(Mk +22.9)2+1. 16 
-0.187(Mk+22.9)2+2.24 
-0.098(Mk+22.9)2+1.09 


-0.223(M,.,o+20.8)2+1.10 
-0.318(Mr,o+20.8)2+2.04 
-0.095(M^,'o+20.9)2+0.94 


1.027(Mk)7(c) 
0.487(Mk)7(c) 
0.947(Mk)7(c) 


1.067(M^,o)7(c) 
0.527(M,,o)7(c) 
1.037(M^,o)7(c) 





* Range of validity: 1.74 < c < 3.06*" Range of validity; -23.25 < Mk < -19.95'^ Range of validity; -21.95 < Mrfi < 
— 18.65* Note: Expressions for 7^ can be obtained by 7^ = ju-K — 7k ~ 7u-K, where we assume jk is small. Earlier 
studies (e.g. TuUy et al. 1998; Masters et al. 2003) showed ^ 0.2. Similarly, we have 77- = 7, — k — Ik ~ Ir-K- 



TABLE 2 
Fitting Functions (for the median) 





Dependence on c 

7(c)- 


Dependence on Mk 


Dependence on M^fi 

l(MrfiY 


Mk kc 
i(cMk) 


Mj.,0 & c 
■y{c,Mr,o) 


u-Ks 

u — r 


-1.13(c-2.49)2+1.08 
-2.84(c-2.61)2+2.49 
-1.70(c-2.69)2+1.43 


-0.072(AfK+23. 1)2+1. 13 
-0.154(A/k+23. 1)2+2.21 
-0.082(Ma'+23.1)2+1.08 


-0.215(Mr.o+20.8)2+1.08 
-0.296(Mr.o+20.9)2+1.99 
-0.081(Mrj)+21. 2)2+0.92 


1.077(Mx)7(c) 
0.527(Mk)7(c) 
1.027(Mk)7(c) 


1.097(Mr.o)7(c) 
0.557(M,.o)7(c) 
1.147(Mr,o)7(c) 





Range of validity: 1.74 < c < 3.06*= Range of validity: -23.25 < Mk < -19.95^= Range of validity: -21.95 < Mrfi < 
— 18.65 Note: Expressions for 7^ can be obtained by 7^ = ■yu-K — iK ~ fu-K, where we assume fK is small. Earlier 
studies (e.g. Tully et al. 1998; Masters et al. 2003) showed 7k ^ 0.2. Similarly, we have 77- = 7, — k — Ik ~ Ir-K- 



